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<D ' Abstract 

Q, 

, We study the stability of some classical string worldsheet solutions employed for computing the potential 

fNj ■ energy between two static fundamental quarks in confining and non-confining gravity duals. We discuss 

the fixing of the diffeomorphism invariance of the string action, its relation with the fluctuation orientation 
and the interpretation of the quark mass substraction worldsheet needed for computing the potential energy 
in smooth (confining) gravity background. We consider various dual gravity backgrounds and show by 
a numerical analysis the existence of instabilities under linear fluctuations for classical string embedding 
solutions having positive length function derivative L'(tq) > 0. Finally we make a brief discussion of 't Hooft 
loops in non-conformal backgrounds. 
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1 Introduction 



The proof of confinement in non-abelian gauge theories from first principles remains to date unsolved. The 
strong-coupling aspect of the phenomenon precludes from attacking it with standard QFT perturbative 
techniques. Nevertheless a criteria for confinement was proposed long ago by Wilson [T|. The criteria states 
that an area law behavior for the so called Wilson loop indicates confinement of the chromoelectric flux 
tubes (QCD string). In particular, an area law behavior for a rectangular (infinite strip) spacetime contour 
corresponding to a static quark-antiquark pair indicates a linear confining potential between quarks. 

In the latest years, the gauge/string correspondence [2] has provided new insights into the confinement 
phenomenon. The crucial observation is that the string aspect of the chromoelectric flux tube manifests 
in the dual perspective by the appearance of a holographic dimension |: 2|] , [3| - A string theory prescription 
for computing Wilson loops was proposed in [4]: the gauge theory loop is to be though at infinity in the 
radial holographic coordinate and the Wilson loop for fundamental quarks is defined by an open string whose 
endpoints lie on the loop at infinity. In the large N c limit ('t Hooft limit) the QCD string self-interactions 
vanish [6] and the dual gravity prescription for computing the potential energy between quarks in a given 
gauge theory amounts to finding a minimal surface in the corresponding gravity dual. 

The canonical computation for the potential energy between a pair of fundamental static quarks (rectan- 
gular loop) involves a U-shaped string extending in the holographic direction where the gauge theory quarks 
separation L translates into the separation between the fixed open string endpoints located at infinity. As the 
quarks separation varies, the string worldsheet explores the holographic direction, therefore, the minimum 
radial position tq reached by the string depend on the endpoints separation distance L. This procedure was 
applied to a number of paradigmatic examples and gave results consistent with gauge theory expectations, 
in particular a theorem stating sufficient conditions for confining backgrounds was proved in [8j (see [7] for 
a review) . Extensions to higher gauge group representations and 't Hooft loops were analyzed and proposed 
in [9J,[10J,[11J, they involve higher dimensional Z?-branes with or without worldvolume gauge fields turned 
on. 

In many applications to dual gravity backgrounds the prescription [4] has been applied at the classical 
zeroth order level to establish confinement, phase transitions or transport properties [7].[12|. only recently 
has the stability of some classical string embeddings been studied [13J,[14J,|15J (see also [17J,[18J,[19|,|20J). 
One of the motivations for the stability analysis, in generalized situations, has been the appearance of 
multiple classical embedding solutions for given boundary conditions [I4j,[15j,[2l] (see also [22]) signaled by 
the presence of extrema in the length function L(ro) (see [23j for related recent work). The presence of a 
maximum separation length was interpreted as dual to the occurrence of screening. The aim of the present 
work is to show that whenever one has a L'(ro) > branch of solutions, they are unstable. We will confirm 
this statement by explicitly showing the existence, in particular gravity backgrounds, of unstable (u 2 < 0) 
modes for the L'(ro) > branches. It is worth mentioning that this is a satisfying result since the expected 
physical behavior for the L(r ) relation from the gauge/string correspondence is to have L'(r ) < 0. We will 
mention briefly an analysis of 't Hooft loops computations in non-conformal gravity duals where instabilities 
also arise. 

Along the way we will discuss various aspects of the classical embeddings: the first one regards the 
physical interpretation of the configuration employed for obtaining a finite potential energy between quarks 
in smooth gravity backgrounds, the second regards the diffeomorphism invariance of the string action and 
its relation to the possible gauge choices for the orientation of the in-plane fluctuations, the third one is 
the relation between instabilities of the string embedding and the L(ro) relation between the separation of 
the string endpoints at infinity L and the maximum depth reached by the string probe in the holographic 
direction r$. It was proved in [15] that the presence of an extremum in the L(r$) relation leads to the 
existence of a zero mode for the longitudinal fluctuations, signaling an instability. We will confirm this fact 
by explicitly computing the lowest fluctuation modes in a number of gravity duals examples. 

The paper is organized as follows: in section [2] we review the prescription for computing Wilson loops 
from gravity backgrounds. In section [3] we describe the backgrounds we will study and compute their length 
and energy functions. In section |4] we perform a quadratic fluctuations analysis and compute numerically 
the lowest fluctuation modes. In section [5] we reanalyze the result of section [4] transforming the fluctuations 
equations of motion into a Schrodinger problem. In section [6] we briefly discuss the 't Hooft loop case and 
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in section [7] we summarize our conclusions. We conclude with two appendices with technical details. 



2 Wilson loops and string solutions 
2.1 Static string U-shaped embeddings 

The starting point for Wilson loop computations, in the large N c , A ('t Hooft limit), from gravity duals with 
metric g^ is the Nambu-Goto actiorQ 



S=- - / dTda\/r]h. (1) 

2ira J 

Here h — deth a 0, h a p — g^daX^d^X" is the induced metric on the string worldsheet, d a = d/d^ a with 
£ a = {r, a} the string worldsheet coordinates and X^ run over the target space coordinates. A sign r\ 
accounts for possible Euclidean (rj = +) and (timelike) Lorentzian (r/ = — ) configurations. The class of 
metrics we consider take the form 

ds 2 = -g t {r)dt 2 + g x {r)dx\ + g r {r)dr 2 + g ab (r, 9)d9 a d6 b . (2) 

The t, Xi (i — 1, 2, 3) coordinates represent the gauge theory coordinates, r is the bulk holographic coordinate 
and 8 a (a, b=l, .., 5) are additional angular coordinates parametrizing a compact 5d space £5. The potential 
energy between quarks involves solving the NG action for (timelike) worldsheets corresponding to strings 
whose endpoints at infinity lie on the loop to be computed, typically the endpoints are kept separated by a 
constant distance L in one of the X4 coordinates which we call x (see [7] for a review) . We start analyzing 
static embeddings of the form t(r),x(a),r(cr), with all other coordinates fixed to constant^]. Placing the 
anzats into the action leads to the correct equations of motion, one therefore has 



S = f dTdaJg t (r)P (g x (r) & + g r (r) r 2 ) 

2na J v 



1 



2na' 

T 
2na' 



dtdayj g t (r) (g x (r) i 2 + g r (r) f 2 



da^p{r)x 2 +g 2 (r)r 2 . (3) 



where g 2 (r) — gt{r)g r (r) and f 2 {r) — gt{r)g x (r). The reparametrization invariance of J!]) factorizes the 
temporal extension of the loop T and reduces the Wilson loop computation to finding a geodesic in the 
effective 2-dimensional geometry 

ds 2 eff = f(r)dx 2 + g 2 (r)dr 2 . (4) 
The conserved charge associated to ^-translations in J3]) is 

_ A 



from which one obtains 

x{a)=±A 9 }j'l . 1 f(a). (6) 

/(r) v // 2 (r) - A 2 

Reparametrization invariance guarantees that solves the r-equation of motion. Calling ro the point given 
by f{ r o) — A, J6]) can rewritten as [8] 

dx =± gix) f(r ) 

dr f(r) v// 2 (r)-/ 2 (r ) • 



1 Generically one should take into account contributions from B2 background fields but for the ansatz we will consider they 
do not contribute to JTJ. 

2 Generalizations considering moving quarks on the boundary, relevant for QGP applications, are straightforward and 
rephrased in terms of non-diagonal terms in the metric[l4].|21|. 
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The boundary conditions at infinity for the string endpoints separation are Ax\ r=00 = L. From J7]) one 
notices that the string reaches the boundary in an orthogonal way. Two natural gauge choices that appear 
in the literature are: x(a) = a (x-gauge) or r(a) — <r (r-gauge). The first choice (x-gauge) has the 
benefit of providing a complete parametrization of the embedding r(x) when imposing x £ [— L/2, L/2] and 
r(±L/2) = oo (the tip of the string is conventionally chosen to be at (r,x) — (r ,0)). Making x — ► t, 
equation (7]) can be understood as a zero energy motion in a potential U(r) given by 



U(r) 



P(r)(f 2 (r)-f 2 (r )) 
S 2 (r)/ 2 M 



(8) 



the point ro is therefore easily seen to be the minimum value in the holographic coordinate reached by 
the string. The second choice (r-gauge) gives a double valued x(r) relation when imposing r £ [ro,oo) and 
x(oo) = ±L/2. Nevertheless in several examples leads to closed analytical expressions for x(r) and moreover, 
when computing fluctuations around the static solution, drastically simplify the equations of motion since 
no fluctuations in the metric components g^ should be taken into account (see [5].[15j). One should keep in 
mind that the tip of the string is a special point in the r-gauge since we must patch there the two branches 
corresponding to the ± signs in ([7]) (see the following section). 

Integrating (J6j) we arrive to the important L(r ) length function, 



L(r ) 



9(r) 



/(ro) 



: dr . 



(9) 



/W VP(r)-f 2 (r ) 

Assuming that /(r), f'(r) > the lower limit ro in ([9} is generically integrabkjf]. Note that a finite lhs in © 
demands that g/ f 2 should decay at infinity faster than 1/r. In the following sections we will be interested 
in the relation ([9]). The derivative L'(r ) can be computed as follows [H] (see also the recent work |33j), 



L'(r ) 



9{r) 



9(r) 



V/ 2 (r)-/ 2 (r ) 
f(r )g(r) 



/'(ro) 
/'(ro) 



f(r)g(r) 



, : (/ 2 W-/ 2 (r ))i 
g(r) d ( 



■ dr 



dv 

fir) dr 



V/ 2 W ~ / 2 (ro) 



f'(rWP(r) - /2(r ) 



dr 



/'(ro) 



V/ 2 (r) - / 2 (ro) dr V/'W 



S(r) 



where we have integrated by parts when passing from the second to the third line. Since the first term in 
the rhs of the third line vanishes in all the backgrounds we consider, we obtain 



L'(r ) = 2 



dr 



/'(ro) 



Vf 2 (r)~P(ro)dr \f'(r) 



g(r) 



(10) 



The energy of the qq configuration was proposed in |4] to be given by the length of the string solution 
([6]) in the effective 2d metric Q, 



E 



1 



2na' 



da^ f 2 (r)x(a) 2 + g 2 (r)r(ay 



The expressions for the energy in the x- and r-gauges, using J6]) are 

1 



E(r ) = 



2na 
1 

ira' 



L/2 dx 
-L/2 /(ro) 



dr- 



g(r)f(r) 



V/ 2 (r)-/ 2 M 



(11) 

(12) 
(13) 



3 A zero /'(ro) = leads to a ro-dependent logarithmic divergence in J9j (see (33] for a generic discussion on divergences in 
the length function). 
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The energy computed by expression l|lip diverges due to the infinite extension of the stringy. The inter- 
pretation for this divergence is that fTTj) contains, in addition to the potential energy between quarks, the 
self energy (mass) of the external quarks [4j. In order to obtain a meaningful quantity and get the poten- 
tial between quarks we should compare (fTTj) with respect to a reference state taking care to substract a r 
independent quantity. It is customary to take the length of a straight string stretching from infinity all the 
way down to the interior of the bulk spacetime along the r coordinate, with all other coordinates set to 
constants, as the "bare" quark mass. Calling r = r m i n the minimum allowed value for the radial coordinate 
in the geometry ((2]), either because of presence of a horizon (e.g. AdS in Poincare coordinates or thermal 
BH backgrounds) or because the spacetime ends in a regular fashion (e.g. Witten AdS soliton, Maldacena 
Nunez and Klebanov-Strassler backgrounds) , the quark mass takes the formal 



2na' 



g(r) dr . 



(14) 



The potential energy between quarks obtained from (fTTj) after substracting the quarks self energy (jT4j) is 



£99 M 



E(ro) 
1 



2 m„ 



9{r)f{r) 



Vf 2 (r)~f 2 (ro) 



: dr 



g(r) dr 



(15) 



Eliminating r Q from ((9)) and lfT5]) we obtain the gauge/string proposal for the potential energy between 
quarks in the planar large 't Hooft limit String (L). In the following sections we will plot this relation in 
several examples and analyze its functional form. For completeness we compute the derivative of lfT5|). one 
has 



KM = 



i 



9{r)f(r) 



Vf 2 (r)-P(r ) 
Using the first line of lfT0|) one obtains |15| 



dr 



f(r)g(r)f(r )f(r Q ) 
(/ 2 (r)-/ 2 (r ))l 



KM 



1 



7 /(n>)£'(ro) 



dE„ 



1 



7 /M, 



(16) 



27ra' v ' v ' dL 2na' 

where ro in the last expression should be understood as the function ro(L) obtained by inverting J9]|. 

We end this section quoting some conditions that must be satisfied by any potential pretending to describe 
the interaction between physical quarks. The so called 'concavity' conditions proved in [24] are 



dV n 

dL >0 > 



d 2 V 
dL 2 



< 0. 



(17) 



This conditions hold independently of the gauge group and the details of the matter sector. The physical 
interpretation of lfl7|) is that the force between the quark-antiquark pair is: (i) always attractive and (ii) 
a non increasing function of their separation distance. From (|9|) . 1(15)1 . 1(16)1 we find that the string proposal 
String (L) gives [22] 



String dE q g dr 



dL 



99 

dro dL 



1 



2na 



7 /M, 



d 2 ^string 

dL 2 



1 



27ra' 



(18) 



The first condition is always met in dual gravity backgrounds since by definition f(r) > 0. Although in 
all our examples f'(r) > 0, the second condition might fail whenever L'(r ) is positive. We will present 
cases where this non-physical behavior appears and show that precisely in those circumstances the string 
embedding solution (Zl)-© is unstable under small perturbations. This last statement was the motivation 
of the present work. 



4 The divergence in Ijl3ll is generically ro -independent, an additional ro-dependent divergence might appear in J9]| and 
whenever the string stretches to regions where f (ro) = 0. 

5 See the discussion at the end of sections 13.1.21 and 13.31 regarding the interpretation of the reference state in smooth gravity 
backgrounds. 
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2.2 Stability analysis of classical string embeddings 

We will study in this section the stability of the classical solution {r c \(<r),x c \(a)) given by I©- ([9]) under small 
(linear) perturbations. A general fluctuation around the embedding solution can be written as 

X» = (r, x d (a) + S Xi (t, a), Sx 2 (t, <t),5x 3 (t, v),r d {a) + Sr(r, a), 9 a + 69 a {r, a)) . (19) 

We can use the diffeomorphism invariance of the action to fix t = t and forget about the t-equation of motion. 
For the class of metrics considered in , the 8x2 and 6x3 fluctuations decouple and as expected satisfy the 
same equation of motion, the S9 a fluctuations mix among themselves for generic compact manifolds and lead 
to five eom (we will not analyze the angular fluctuations in the present work and we consistently set them 
to zero), finally the 5x 1 and Sr fluctuations result mixed in two coupled equations. 
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Figure 1: r-gauge fixing: The dashed line rep- 
resents the classical embedding over which we 
perturb. At the tip ro, the fluctuation is ori- 
ented along the string, and therefore not phys- 
ical. 



Figure 2: a;-gauge fixing: The dashed line rep- 
resents the classical embedding we are per- 
turbing. This is a physical gauge choice 
all over the embedding solution, the plotted 
(even) fluctuation changes the position of the 
tip. 



It is easily shown that the r- and xi-equations of motion are proportional. The remaining diffeomorphism 
should therefore be used to fix the orientation of the (Sr, 5x\) vector (in-plane fluctuation) at each point 
of the solution (r c i(er), x c \(a)). After imposing a gauge constraint one equation describes the fluctuations 
in the (r, xi)-plane and we end up with a well posed system of differential equations. The physical gauge 
choice (n-gauge) would be to orient the fluctuation along the normal direction to the classical embedding 
(r c | (<r), x c \ (a)) and generically this means to account for fluctuations on both x\ and r coordinates. Two other 
natural possibilities considered in the literature correspond to fixing Sx\(t, a) — (x-gauge) or Sr(r, a) = 
(r-gauge). The n- and x-gauge fixings, as mentioned in the previous section, parametrize the fluctuations 
along the whole classical embedding but the equations of motion are lengthier because the Sr fluctuations 
result in additional contributions to the eom arising from changes in the metric |18J.[15J. Note that at first 
sight the n- and x- gauge fixings appear to allow for the oscillation of the tip, while the r-gauge fixing should 
not (more on this in the following). We will choose to work in the r-gauge, this means fixing 5r(r,a) = 
and work with simpler equations defined on half of the embedding. It will then be mandatory to analyze 
the boundary conditions to be imposed at the tip (ro,0) of the embedding in order to get a meaningful 
solution. Moreover, since we are considering fluctuations along the x\ coordinate, precisely at the tip the 
5x\ fluctuation is oriented along the string worldvolume and not transverse to it: the 5x\ fluctuation in the 
r-gauge is therefore not physical at the tip of the embedding. The required additional analysis developed in 
[14], [15] will be discussed below. At last, another advantage of the r-gauge is that it gives closed expressions 
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for the linearized fluctuations equations of motion (see l(2"2"|) - ([2"3"10 . the isometry along X\ implies that upon 
computing fluctuations over the classical solution (7|)-((9|), we do not need the explicit analytic solution x c \(r), 
what contribute to the equations of motion is its derivative (O (as an example compare (|113p with (|114p ). 
We would like to quote the work [14] where Sr fluctuations in the 5x\ = gauge fixing were considered over 
(half) the classical embedding parametrized in the r-gauge (see below eqn. ([27 ]) . ([28]) ). 

In the following we will be concerned with the equations of motion for the 5x{ fluctuations, they can be 
obtained by placing the ansatz 



t = t, xi = x c ](r) + 8xx(t, r), x 2 = Sx 2 (t,r), X3 = Sx3(t,r), 
into the action (Q]). Expanding to second order in the fluctuations one obtains 



(20) 



1 



h 2 (r) (f(r) - / 2 M) (Sx 1 ) 2 - (f(r) - f 2 (r )) 2 (Sx r) 5 



g(r)f(r)VP(r)-f 2 (r ) 
+f(r)h 2 (r)((Sx 2 ) 2 + (Sxs) 2 ) f(r) (f(r) - f 2 (r )) ((Sx 2 ) 2 + (5x 3 ) 2 ) 



(21) 



where h 2 (r) = g x (r)g r (r). The Euler-Lagrange equation for the 5x\ fluctuation is 



(P(r) - f 2 (r ))i d 



g(r)f(r) 



dr 



2 h 2 (r)^f 2 (r)-p(r ) 



9{r)f{r) 



where we factorized the time dependence of the fluctuations as Sx(t, r) 
fluctuations transverse to the (r, x\ )-plane obtained from ([21]) are 



Sx r (r)=0, (22) 
6x(r) e~ lujt . The equations for the 



/(rV/ 2 M-/ 2 M d 



h 2 (r)f(r) 



9(r) 



dr 



Sx m (r) = , m = 2, 3 . 



(23) 



Equations |22|) - (|23| are differential equations of the Sturm-Liouville type defined for the half line r min < 
r < r < oo and we are interested in analyzing the existence of instabilities, in particular determining the 
range of values of r$ for which u> 2 < 0. 

The boundary conditions to be imposed on the problem are Dirichlet, this means fluctuations keeping 
the string endpoints fixed at the boundary Sxi(r, a)\ r=QO — 0, but, the nature of the r-gauge parametrizing 
only half of the (r c \(a), x c \(a)) embedding requires an additional analysis of boundary conditions at the tip 
r = ro (singular point of l(^2 |) - (|23|l ). We start by analyzing ([23|) . the expansion around the tip givefl 



d_ 

dr 



d5x m (r) 



dr 



Lo 2 h 2 {r ) 



2/(m)/'M Vr~^ 



Sx m (r)^0 =S> Sx m (r) w C + C 1 y/r - r + 0{r - r ) . 

(24) 

Here Cq,i are arbitrary constants corresponding to the two independent solutions of the differential equa- 
tion (|23p. which once chosen determine the whole series expansion for Sx m (r). Physically they correspond 
respectively to even and odd fluctuations around the tip once we patch them with the fluctuations around 
the other half of the embedding, which obviously satisfy the same eom. A discrete set of eigenvalues uo 2 is 
expected if non-normalizable solutions exist in the large r-limit, we also expect the even solution to have the 
lowest uj 2 eigenvalue (see appendix IB} . 

We now turn to the analysis of the in-plane Sxi fluctuation. Expanding (|22| around r = rg one finds 



d_ 

dr 



(r - r ) 



3 dSxi(r) 

2 L-L 

dr 



u} h (ro) 
2/(r„)/'(r 



■y/r — r Sxi(r) w 



8 Xl (r) e*C' Q + C[ 



1 



\/r - r 



+ 0(Vr~Vo-). 

(25) 

A singular behavior appears for 6x% at the tip and one might be tempted to cancel it by imposing C[ = 0. 
Nevertheless, as mentioned above one should take notice that the r-gauge fixing implies that, at the tip, Sxi 
is directed along the string worldsheet and not transverse to it, therefore in the r-gauge the displacement 



3 We consider f(r) to be an increasing function of r having no zeroes except perhaps at the bottom of the bulk r : 



Sxi at the tip is not physical. In order to give a physical interpretation to the constants C 01 in ((25j) we now 
make a change from the r- to the x-gauge (cf. [12], [15]). The change of gauge on the ansatz 1|20[1 can be 
implemented by a change of variables on the solution (|22|) from r to a new variable which we call u. It can 
be implemented perturbatively, to first order in the fluctuation the relation is [15] . 

u = r + A(t,r) with Afar) = 5xi ^ , (26) 

x' d (r) 

This transformation performs the desired change of gauge since 



%i = x d{ r ) + Sxi(t, r) — x c \(u — A(t, r)) + 8xi(t, r) 
u) - x' d (r 
Sxi(t, u) 



x d (u) - x' d (r) Sxi ^'^ +6xi(t,r) = x c \ (u) (27) 



x' ci (u) 



(28) 



here r$ < u < oo. This is precisely the gauge fixing employed in [14] mentioned above. The second term 
on the rhs of (|28|) is interpreted as the r-direction fluctuation induced by the (r-gauge) xi-fluctuation. It 
is now easy to see that it is finite. The asymptotic behavior lj25]) and the expansion around the tip of {7]), 
x' d (r) ~ (r — r ) -1 / 2 , plugged in (|28| gives 



r » r - a(C' V^^ + C[) + 0(u - r ) , (29) 

with a a finite constant. The result (|29| shows that the physical Sr fluctuation originated from the non- 
physical divergent Sx% fluctuation at the tip is manifestly finite. Therefore we interpret the (Cq) C[ in (|25|) 
as the fluctuation that (do not) oscillates the tip position. 

In the following sections we will show numerical solutions of |22|) for various gravity backgrounds deter- 
mining the lowest eigenvalues leading to normalizable solutions. We will solve (|22"1) by a shooting method 
integrating numerically from the tip tq up to a large value Too. The allowed values for ui 2 will be obtained 
imposing the numerical solution to be zero at r = r^. The boundary conditions at the tip corresponding 
to even solutions are Cq — and C[ arbitrary, for numerical purposes we set C[ = 1, its value sets the 
normalization of the fluctuation. An even solution around the tip satisfies 



dSr(t, r) 



dx\ 



= where Sr(t,r) = _ fe H*' r ) , ( 30 ) 

x d( r ) 



Using (|27| we can write f30|) in terms of 6xi(r). The boundary conditions for even solutions of (1221) are 
implemented numerically as 

5x 1 (r) + 2(r-r )^rl=0, r - r 
ar 



y/r- r Sxi(r) = 1, r -> r . (31) 
Odd solutions C = 1 and C{ = are implemented as 

bxi{r) + 2(r - r ) — = 1, r -> r 



Vf^fei(r) = 0, r^r . (32) 

Summarizing, in general backgrounds, the functional relation of the classical solution ([7]) between the x\ and 
r coordinates at the tip takes the form x d (r) w r — r Q , and the asymptotic behavior of the Si-fluctuations 
over it, in the r-gauge, is of the form (J25j) . Although a divergent piece appears in (|25|) . an appropriate change 
of gauge shows that the divergent and non-divergent pieces correspond respectively to (physical) even and 
odd fluctuations around the tip. 
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3 Gravity Backgrounds 



In this section we compute the string embeddings ([6]) dual to rectangular Wilson loops for a a number of 
paradigmatic gravity backgrounds. We review the AdS§ x S 5 [4j and AdS's-SchwarzschildxS' 5 [16] cases. Next 
we perform the numerical analysis of the equations iJH) and lfl5"Tl for the Maldacena-Nunez [25] , Klebanov- 
Strassler [26] and the generalized Maldacena-Nunez [28j backgrounds. In all cases the geometry is supported 
by some non-trivial p-form fluxes, but they will not be relevant for our computations. 



3.1 AdS 5 x S 5 

This background is dual to J\f = 4 SYM with G = SU(N) gauge group in the Coulomb phase. The 
AdS curvature R relates to the gauge theory 't Hooft coupling A as R A — a' 2 A and the flux of the 5- 
form supporting the geometry N — f s5 F 5 relates to the rank of the gauge group as TV = Rank(G) [2]. 
The conformal invariance of the gauge theory implies a Coulomb behavior for the potential V(L) ~ 1/L 
between quarks. The novelty of the gravity computation is to compute the gauge coupling dependence of 
the proportionality coefficient. 

3.1.1 Poincare coordinates |4j 

This coordinate system is supposed to describe the gauge theory formulated on R 3 ' 1 . The metric is written 
as 

ds 2 = -^{-dt 2 + d Xl dxi) + R 2 — + R 2 dn 2 r . (33) 
R z H 

One finds f 2 {r) — r 4 /R 4 , g 2 {r) — 1. The radial coordinate range isO<r<oo, atr = one finds a Killing 
horizon. Equation J7]) can be analytically solved in the r-gauge, one obtains 



-c,W=±{c*e-^B((^) 4 ;|0} ) , <,<oo 



(34) 



where B(z; a, b) is the incomplete beta function B(z; a, b) = J* t a 1 (1 — t) b x dt. The boundary conditions fix 
the constant in f34|) and relate the parameters ro and L, setting x c \(ro) = and ar c |(oo) = ±L/2 one obtains 

®, 

R 2 n /3 1\ R 2 (2tt)! 
^V'o; = ■ - 

The energy (f]~5jl takes the form [4] 



^»"5; B U'5r^W- (35> 



, 3 



^^^^H-M-^S- 06) 



L 4 J 



here K(m), E{m) are the complete elliptic integrals of first and second kind. Eliminating ro from expressions 
(|35| - l|36p the AdS/CFT proposal for the interaction potential between fundamental quarks in the large 't 
Hooft coupling for the TV = 4 SYM theory is [4] 

(2n) 2 R 2 /a> VA 
Vstnng(-£<) — p^Tjl f, T ■ (3?) 

An attractive Coulomb potential is obtained as expected from conformal invariance. The interesting result 
is the VA = {9ymN)2 dependence when compared with the perturbative A = gy M N result. This suggests 
that some renormalization of the charges takes place at strong coupling [4J. Note that in order to obtain 
the negative quantity l(37|) starting from a positive definite one (see eqn l[TTj) ). it was crucial to substract the 
quark masses lfT4|) . 
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3.1.2 Global coordinates 



We discuss this example because it clarifies conceptual issues regarding the interpretation of the substraction 
procedure lfl5l) in smooth and complete gravity backgrounds (see sections 13.31 13.41 13. 5ft . 

Computations performed in AdS global coordinates are supposed to represent the Af = 4 SYM gauge 



theory defined on S 



The AdS metric is now written as 



ds 2 = R 2 [~ cosh 2 pdt 2 + dp 2 



sinh 2 pdVL 2 ] 



(38) 



All coordinates are dimensionless in this case with the AdS radius R setting the scale. We write the S 3 
metric as dQ 2 = d9f + sin 2 0i(d#2 + sin 02dtp 2 ). Being tp a cyclic coordinate, the appropriate string ansatz 
is t = r, p = p(a), tp = <7, one obtains f 2 (p) = \ sinh 2 2p, g 2 (p) — cosh 2 p. The remaining angular variables 
must be set to Oi = ? (equator of S 3 ) in order to satisfy their eom. The conserved charge in the ip coordinate 
leads to an effective one dimensional zero energy motion 



p 2 + U{p) = , 



(39) 



where p — dp /dip and the potential U(p) = sinh 2 /? ^1 — ^^r§^j ■ Po is the minimum radial position reached 

by the string when separating the string endpoints at infinity by Atp = $. The $(po) relation © is computed 
straightforwardly 

/ \ f°° sinh2p , , N 

*(A>) = 2/ w , =dp. (40) 

Jr sinhpv sinh 2p — sinh 2po 

Since the gauge theory is defined on a S 3 there exists a maximum separation for the quarks and it corresponds 
to placing them at antipodes on the equator of the S 3 . This results in the string reaching the origin $(0) = n 
and leading to a smooth straight worldsheet (stretched along the radial r coordinate) parametrized by two 
halves at tp — <po and tp = it + po . The (divergent) energy lfl3l) of the configuration lj39j) is 



E(po) 



R 



sinh 2p 



2na ' Jpo sinh py/smh 2 2p - sinh 2 2 / 9 



dp . 



(41) 



Substracting the quark masses as in lfl5|) leads to, 



R 
2na' 




2 cosh p 



sinh 2 2pn 
sinh 2 2p 



2 cosh p\ dp — 2 sinh po 



(42) 



which is finite and negative definite (see figE]). The finite result |42|) should be understood as resulting from 
the comparison of l(4T|) wrt the aforementioned smooth reference state consisting in a straight string with 
its endpoints at infinity at antipodes on the S 3 equator. We interpret this last configuration as the one 
corresponding to "infinitely" separated quarks on S 3 . Note that the reference state and the configuration we 
are analyzing satisfy different boundary conditions. 



3.2 AdS'5-SchwarzschildxS' 5 

Finite temperature gauge theories are described by considering black hole (BH) solutions in the gravity duals 
[5]. The near horizon geometry of N non extremal (black) D3-branes is therefore conjectured to describe 
N = 4 SYM at finite temperature. As explained in [5], the appropriate BH background for describing the 
gauge theory on R 3 ' 1 involves a delicate infinite mass limit of the ^4dS5-Schwarzschild BH resulting LqJ, 



ds 2 = — 
R 2 



d-^)dt 2 

r 4 



dx 7 dx I 



R 2 



1 



-dr 2 



R 2 dn 2 . 



(43) 



7 The limit results in the metric depending on only one scale, contrary to the finite mass case where one has two parameters: 
the temperature (BH mass) and the AdS radius. This last geometry was shown to exhibit a phase transition (Hawking-Page), 
which was interpreted as dual to the confiniement/deconfinement phase transition in Af = 4 on 5 3 [5]. 
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Figure 3: String obtained from f40|) and (|42|) . For small angular separations 
between the string endpoints $ -C 1 one finds the expected Coulomb behavior 
V ~ \/A/(.R$). For larger separations the solution deviates due to the compact- 
ness of the S 3 . 



The BH horizon is located at r = fi and its temperature is T = ^s-7r. It is convenient to work with 
dimensionless coordinates, scaling r = fx p, t — R 2 /pt and x = R 2 /py one obtains 



The scale of the dimensionless gauge theory coordinates t, y, in (|44|) is set by i? 2 /V and one finds f 2 (p) = 
p 4 — 1, g 2 (p) = 1 and p = 1 as the horizon location. The expressions for the dimensionless qq separation 
length ([9]) and potential energy lfT5|) can be analytically computed jlol \W\ [22] 



Here po > 1 is the minimum radial position reached by the string and the minimum radial value r m i n in 
(fl5l) was taken to be the horizon location r min = p. One can easily check that in the small temperature 
limit LT <C 1 (corresponding to po 1) one recovers the zero temperature behavior (J35J) - (J3TJ) . We have 
plotted in figures |4] and [6] the behavior of the length (|45l) and the energy f46|) as functions of po- In figure 
[5] we plotted the relation V 5 t rmg (L) obtained from f45|) - (|46| by eliminating po, the result is a double valued 
function. 

Figure [4] shows a maximum L c ~ 0.869 at po c — 1.177 which implies that no smooth solution connecting 
the pair of quarks exists for L > L c . The only existing solution for L > L c corresponds to two straight 
strings reaching the horizon. This configuration, used for the substraction in lfT5|) . is interpreted as the one 
corresponding to a pair of free quarks. The existence of a maximum in the L(ro) relation in BH backgrounds 
has been interpreted as the gravity dual of thermal bath screening [16]. Figure [4] also shows the existence 
of two branches of solutions for each string endpoints separation distance L < L c . The left branch [L 1 > 0) 
leads to a potential V^ tr i ng (L) not satisfying the conditions lfl7|) and it should therefore be non-physical (upper 
curve in fig. [5]). The rewarding result as we shall see is that string theory is wise: the left branch should not 
be trusted since it is unstable under small perturbations (see section [472)) . 

The last point to comment is that although, as we shall show, all the L' < region (1.177 < po < oo) is 
stable under small perturbations, one expects the lower curve in figure [5] solution to be metastable whenever 
E qq > (0.754 < L < 0.869 or 1.177 < po < 1.524). The reason for this is that the (reference) two straight 
lines solution has E q q — and is therefore the absolute stable minimum fro 0.754 < L < 0.869 (see fig. [6|). 
Taking this last fact into account [IB] suggested to take L max ~ 0.754 as the screening length. 




(44) 




(45) 




(46) 
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Figure 4: L(po) relation (|45|) : A maximum is ob- 
served for p — 1-177. The two branches to the 
left /right of the maximum at po c ~ 1.177 lead to 
a double valued V s t r \ ng (L) relation in fig. (5) 




Figure 5: Double valued V^ tr i ng (L) obtained 
from l(45|) and 1(46)1 eliminating numerically 
Pq. The upper curve corresponding to the left 
branch in fig. |4]does not satisfy the conditions 

G3- 
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Figure 6: E q g(p ) relation for thermal Af = 4 SYM. Eqn. 1(16)1 guarantees an extremum in E q q(r ) 
and L(r ) at the same value of r . The curve cuts the axis at po m ~ 1.524, for p > p 0m the U- 
shaped solution is an absolute minimum, for po < po m the two straight strings (reference) solution 
is the absolute minimum. 



3.3 Maldacena-Nunez background 

The r ss of this background is supposed to describe qualitatively the IR regime of d = 4 Af = 1 SYM 
theory. The probe brane configuration leading to this solution consists of TV Z?5-branes wrapped on a finite 
2-cycle at the origin of the resolved conifold. When the backreaction of the D5 on the geometry is taken 
into account, a transition flop occurs (see [30]) leading to a geometry with a smoothly collapsing S 2 and a 
finite S 3 at the origin as for the deformed conifold (see [31] for a review of the solution) . The solution was 
independently found in the context of gauged supergravity in [32j. 
The metric can be written as 1 25 1 



ds 2 = a'Ne 4 



dt 2 + dxidxi + dr 



2 ' e 2h (d9 2 



+ sin 2 Qd<p 2 ) + ^(w l -A 1 ) 2 



where w l (i = 1, 2, 3) are the su(2) right-invariant forms 

w 1 +iw 2 = e^ i ^(d9 + ism0d(p), w 3 = d^ + cos8dip, 



(47) 



(48) 
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and A 1 , 4>, h are given by 
A 1 

o(r) 



,,2/i 



,20 



-a(r)d6, 

2r 
sinh 2r 

rcotli2r — 



A = 



sin 6 dip , A 3 = — cos 9 dip 



sinh 2r 



sinh 2r 

2e' 1 



(49) 



4>o is an integration constant which sets the value of the dilaton <j> at the origin (some authors write g s = ) . 
The t,Xi,r coordinates are dimensionless and its scale is set by (a'N)?. The metric Q47jl is regular in the 
r — > limit with the 2-sphere (6, <p) smoothly collapsing and the resulting topology of the spacetime is of the 
form M7 x S 3 , contrary to the cases discussed before where the spacetimes were of the form M5 x X 5 with 
X 5 compact. The reason for this is that the background lf47|) models the 5+1 gauge theory on the wrapped 
Z?5-branes. Nevertheless, one expects to get an effective 3+1 theory at energies E < l/Rgi, where i?g2 is 
the radius of the sphere wrapped by the D5. 

The paradigmatic static U-shaped ansatz leads to / 2 (r) = g 2 {r) = e 2 ^ and the separation length |(9|) and 
energy 1|15|1 take the form 



L(r ) 



OO 



E, 



( \ N 

7T 



=20(r o ) 



: dr 



: dr 



, ro Ve 2 ^ - e 2 ^ r «i Jo 
The expression for the energy (|5l|l can be rewritten in an illuminating form \8\ 

. e U(ro) _ e 20(r o ) 



(50) 
(51) 



E qq( r 0) 



N 


[f( 


7T 




N 


e <P(ro) 1 




7T 





Mr) 



dr 



^ r) di 



^ (r) dr 



(52) 



In the large L limit the string reaches the bottom of the bulk (tq — > 0, see figure [7]) and the last two terms 
in lj52j) do not contribute (see [8]). In the large L>1 limit we obtain reinserting units [25] 



String (£) 



27ra' 



-L 



-^string — 



2?ra' 



(53) 



The background (|47|) is therefore predicting linear confinement for large quarks separation in agreement with 
its proposed d = 4 TV = 1 SYM dual gauge theory. From the (chromoelectric) string tension f53|) we see that 
the value of the dilaton at the origin 0o relates to the dynamically generated scale of the dual gauge theory. 

The divergence L — > 00 at ro — > 



We plot in figures [7] and [8] the £(ro) and £(ro) relations (f50|l and |(5T 
in fig. [7] is due to ^:( e2 ^ <r '')| r ._Q — ( see [SSI f° r a recent discussion). 

The result for the potential V s tr\ns(L) in figure [9] is rewarding, but the linear behavior occurs for configu- 
rations having energies above zero. A concern arises as to whether we should trust the result for V^ring > 
(see last paragraph of section l3~2|) . The substraction in (|5l|l corresponds to a pair of straight strings running 
along the radial direction with the remaining spatial coordinates kept fixed. Being the background regular, 
the strings cannot end at any point in the interior and the only possibility for a smooth reference solution 
is to place the string endpoints at antipodes on the ip coordinate and having both the same x\ coordinate 
(see the horizontal blue line in fig. [inland the related discussion in sect. 13. 1 .2|) . We conclude that although 
the linear confinement occurs for configurations having E qq - > 0, the solution is stable and it cannot decay 
to the reference E q q = state since the two configurations being compared in ([5Tjl satisfy different bound- 
ary conditions. Figure [10] depicts the relevant worldsheets for the rectangular Wilson loop computations in 
smooth gravity backgrounds. 
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Figure 7: L(r ) relation ((50]). MN solution. 



Figure 8: E(r ) relation JH])- MN solution. 



V: 



String 
lr 




Figure 9: Potential energy V str - m g{L) for the qq pair obtained for the Maldacena-Nunez solution 
eliminating numerically ro from l(50|) - ([5T|) . Note the change of behavior from Coulomb like to 
linear. 



3.4 Klebanov-Strassler background |26j 

This background describes a N = 1 quiver gauge theory with bifundamental matter fields transforming 
under SU(N + M) x SU{N). The probe branes configuration leading to this geometry is constructed as: 
N D3-branes in the apex of the singular conifold plus M -D5-branes wrapped in the topological S 2 of the 
conifold and sharing the remaining three dimensions with the D3. The solution is supported by a constant 
dilaton <fi = 4>o > which can be set to have j s < 1 everywhere (contrary to the MN solution) . The Af = 1 flow 
SU(N + M) x SU{N) — » SU(M) realized through a cascade of Seiberg dualities in the gauge theory (see 
[34]) manifests in the geometry in a varying 5-form flux. 
The metric reads [26] 



ds 2 = g s a'M [h ^(r)(-dt 2 + dx^dxi) + h$(r)ds$\ 



(54) 



The deformed conifold metric dse can be written 





' (rfr 2 + (g 5 ) 2 ) 
3K 3 (r) 



+ cosh 2 L ((ff 3 )2 + (3 4 )2) + sinh 2 ((5 i )2 + (g2)2) 



(55) 



where 



K(r) 



[sinh(2r) - 2r]s 
2 3 sinh r 



(56) 
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Figure 10: Geodesies employed for the computation of rectangular Wilson loops in smooth gravity back- 
grounds. The black curved vertical line depicts the U-shaped solution, the horizontal blue line is the smooth 
reference state with respect to which we compare the energy of the black line. The black and blue configu- 
rations satisfy different boundary conditions. 



and the <?i defined by 



e 1 + e 3 



V2 



V2 ' 



v/2 ' 



e 2 +e 4 
V2 ' 



9 5 = e\ 



with 



(57) 



- sin 0\ d(pi, 



d0i 



sin ip d02 + cos ip sin 02 d<p2 , 



e 4 = cos tp d02 + sin tp sin 02 d<p2 , e 5 = dip + cos 0\ depi + cos 02 d(p2 ■ 

The coordinates in lj54j) are dimesionless, the gauge theory coordinates t, Xi scale is set by 
scale of the holographic coordinate r is set by £ c f- The h(r) factor takes the form 



g s Mu 

let 



h(r) = 2* 



dx- 



: coth : 



sinh 2 x 



(sinh2a; — 2x)" 



(58) 
and the 

(59) 



The background (f54|) is supported in part by a non-trivial but the embedding we are considering gets no 



contribution from it [27]. The functions in Q are given by f 2 (r) = 
expressions for the length ([9]) and the energy (fT5| are 

dr h(r) 



l 

h(r) 



E q q(r ) 



L(r ) = 2 
9sM 



V6K(r) ^h(r ) - h(r) 



g 2 {r) = 6K l( r \ ■ The dimensionless 



(60) 



dr 



y/Wo) 



dr 



(61) 



V6K(r) y/h(r ) - h{r) Jo V&K(r) 
In figures [11] and [12] we have plotted these two last expressions. As in the MN case a divergence is expected 
at ro = since dh j^ = Eliminating numerically tq from (|60"ll - l|6ip we plot in figure [13] the V strmg (L) 

r—0 

function. A linear relation for the interaction potential is observed for L ^> 1. Proceeding as in lj52j) we find 
the confining string tension to be 

i e 



string 



2wa' g s a'M \fho 



(62) 
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Figure 11: L(ro) relation <f60|l . KS solution. 



Figure 12: E gg -(ro) relation l[61j), KS solution. 



where ho = h(0) ~ 1.1398. As in the Maldacena-Nunez case, the dominant contribution to the minimal 
area ([J), in the large L limit, comes from the r w region. Again, a concern arises regarding whether one 
should trust the E q q > configurations, as discussed at the end of the last section and pictured in fig. [TO] 
the E q q > embeddings are classically stable. 



3.5 Generalized Maldacena-Nunez solutions |28j,|29j 

This class of backgrounds was obtained in [28] generalizing the solution described in section 13.31 The 
solutions were thoroughly discussed in [29] an interpreted as dual to minimally supersymmetric gauge theories 
containing irrelevant dimension six operators. The operator drastically changes the UV behavior of the 
theories taking the solution 'away' from the near horizon of the Z?5-branes which generate the geometry. 
The analysis in [29] shows that the solutions asymptote, for large r, four dimensional Minkowski times the 
deformed conifold. 




Figure 13: V s t r \ ng (L) relation for a rectangular Wilson loop in the KS solution. 



The generalized MN metric reads 



ds 2 




(u>3 + cos 9dtp) 2 



(63) 
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Figure 14: L(po) relation (|68|) . General- 
ized MN background with blue, violet, yel- 
low and green curves corresponding to (i = 
— 1.8, —1.5, —1, —.8 (poo = 7). A minimum 
quark separation length is observed. 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

Figure 15: E(r ) relation ((69j) . gMN back- 
ground colors as in fig. HH 



Making the change of coordinates dp = e k ^ r 'dr, the functions a, k and / satisfy 



d p a = 



d p k = 



-2 



1 + 2pcoth2p 



2k (acosh2p - l) 2 . 
e r-r+ — + a (2p - a smh 2p) 



sinh 2p 



2(1 + a 2 - 2acosh2p)" 



-1 + 2pcoth2p 



e a sinh 2p (a cosh 2/? — 1) + (2p — 4a p cosh 2p + — sinh 4p) 



d P f = - 



1 



(1 - acosh2p) 2 (-4p + sinh4p) 
(1 + a 2 - 2acosh2p)(-l + 2pcoth2p) 



4 sinh 2p 

and the g{p), h{p) functions in lj63j) are given by 

e 2g 



ocosh2p — 1 2h 
acosh2p — 1 ' 



(2a cosh 2p - 1 - a 2 ), with b(p) 



2p 



sinh 2p 



(64) 



(65) 



The first two differential equations in f64|) have a one parameter family of regular solutions. For small r one 
finds [28] 

4 20 + 36p + 9p 2 



a(p) = 1 + pp + 



a 2fc(p) 



6 + 3/i 



15(2 + a*) 



-P 



(66) 



with p taking values in the interval (— 2, — |). Inserting l|66p into the third equation of (|64| and into ([65 
one obtains 



6 + 3p 



4p 2 



6 + 3p 



2/(p) 



(2 + p) 2 2 



(67) 



The arbitrary constant for / following from (|64|) was factored out as g s in (j63|) . The limit values for p give 
known solutions: the p = — | case reproduces the MN solution of section I5~3l with 4> = 4/ (fc = const.), and 
the case p = —2 case leads to 4-dimensional Minkowski space times the deformed conifold (</> being constant 
in this case). Finally, the p — > oo limit of all solutions (except p, = — |) asymptotes the deformed conifold 
metric (see [28] for details). 

The length |(9|) and energy l[T5]l expressions are given by 



L(po) = 2 



=,4/(po) 



y/ e Sf( P ) _ e 8/( P0 ) 



s fc(p) 



f/p 



(68) 
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Figure 16: Double valued V s t r i ng (L) relation for 
gMN background with p = —1. The upper red 
curve is the unphysical branch corresponding 
to string embeddings to the right of the min- 
imum in figure HH The curve does not reach 
the origin, manifesting a minimum quark sep- 
aration length. 




Figure 17: Normalized linear V^tring(^) rela- 
tion for the physical (left) branches of figure 
fT4l Colors as in fig. HH 



c 8/(p) 



(69) 



Note that the radial integration in both expressions extends up to a finite distance poo. The reason, noted 
in [28], being that ([68)1 is UV divergent We should therefore have in mind that the string computation 
corresponds to probing the dual gauge theory with very massive for p^ 1 (but not infinite mass) quarks. 
Moreover since the string endpoints are fixed at a finite radial distance it can be checked at the probe string 
does not reach the gauge theory brane along a normal direction. 

In figures [14l and [T5l we plot numerical solutions of 1|68|1 and (|69|) for various values of p.. Figure (fT4| shows 
that L(po) attains a global minimum for all values of p (except p ^ — §). In other words, no solution exists 
for quarks separations L < L m i n . It is interesting to note that the minimum value for L is attained, for all 
values of p, in a rather small region in the p coordinate near the origin. Based on the concavity considerations 
discussed at the end of section 12.11 we expect the string worldsheets to the left of the minimum p c to be 
physically meaningful (stable) and the ones to the right of p c to be unphysical (unstable). We will show in 
the following section that this is in fact the case by analyzing quadratic fluctuations around the solutions. 
Negative eigenvalues appear for the L'(po) > branch of solutions. 

We finally plot in figures (fIB"| and (fl7|) the String (L) relation for the generalized MN backgrounds where 
a linear confining behavior is observed. Figure lfl6|) shows the double valued Ktring(i) relation for p, = — 1, 
the upper red branch (unphysical) corresponds to the string configurations to the right of the minimum in 
fig. I|14p . In figure [17l we show the V sirmg (L) relation for the physical branches of figure [14l for several values 
of p. Proceeding as in f52|) one finds that all solutions lead to a /i-independent string tension 

T 5trmg = ^. (70) 



4 Stability Analysis 

In this section we study, for the backgrounds presented in the previous section, the eigenvalue problem given 
by the equation of motion (|22|) for in-plane fluctuations in the r-gauge. We are interested in searching 
for unstable modes. For the reasons discussed at the end of section [2~T| our aim is to show that negative 

8 See |33j for a discussion on divergences when computing the length function L(tq). 
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eigenvalues (uj 2 < 0) appear for string configurations belonging to regions where L'(ro) > 0. We study even 
solutions, this means that we choose C[ = 1 and C — in (|25| as the initial condition at the tipQ We 
numerically implement this conditions as (see l[3"T])) 



fax r +2 r-r , n y = 0, 



V?" - fo fei(r) = 1 , r — > r . (71) 

Solving numerically, the allowed eigenvalues oj 2 for ([22]) are obtained by demanding (5xi(r) to be a normal- 
izable solution 

£xi(r)=0, r^oo. (72) 

For completeness we recall now the relation between zero modes for in-plane modes and critical points of 
the L(ro) function |15j . The zero mode solution of ([22]) can be immediately written down 



5x ( ° } (r)=C 



dr- 



g(r)f{r) 



(f 2 (r) - / 2 (r ))§ 



c 



(73) 



where C, C are integration constants, C — to get a normalizable solution and we set C — 1. Integrating 
by parts in ([73]) and using ifTO]) one obtains 



fa< 0) (r) = - 



dr 



9(f) d I 1 
/'(f) df l vv // 2 (f)-/2(r ) / 
ff (r) L'(r) 



f'(rWP(r)-p(r ) 2/'(r) 
Expanding this last expression around the tip r — one has 

(0),_^_ sM 1 , i'(r ) 



fop (r) 



+ 



V2(/'(r ))l 2/'(r ) 



(74) 



(75) 



Generically the first factor in the rhs of ([75]l is non-zero, so a necessary and sufficient condition for obtaining 
an even zero mode solution (see (71])) requires the second term in ([75]) to cancel, equivalently ro must be a 
critical point of the L(ro) length function [15] . 



4.1 AdS 5 

The in-plane fluctuations equation of motion ([22]) for the AdS spacetime in Poincare coordinates ([33]) takes 
the form 

d I (r 4 — r 4 )* d 



dr 



dr 



8x\{r) = < ro < r < 00 . 



(76) 



Dilatation invariance implies that one should be able to factor out the r dependence. Making r — r p one 
obtains 

V-i)*d 



d 

dp \ p 2 dp I 
The asymptotic behavior (p — > 00) of ([77]) reads 



Sxx{p) =0. 



dp ^ dp 



2 D 4n 



UJ Z R 



whose solutions are 



fei(p) w a H — o- , /)>1. 



(77) 

(78) 
(79) 



3 Even solutions correspond to arbitrary C[ at the tip, its value fixes the normalization of the solution. 
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Figure 18: Lowest numerical eigenvalue uj 2 of 
([76| giving a normalizable solution as a func- 
tion of r . 
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Figure 19: Properly normalized lowest eigen- 
values of fig. [18] as a function of ro (see (|77|) ). 



with a , ai integration constants. The behavior 1(79]) implies that normalizable solutions (a = 0) will be 
found only for particular (discrete) values uj 2 ^. 

As a test for our shooting method, we have numerically integrated l|76p for different values of ro and 
determined the minimum to 2 eigenvalues leading to a normalizable solution. In figure [18] we plot these uj 2 
as a function of ro. They are positive for all values of ro, signaling the stability of the U-shaped string 
configuration. In fig. [19] we show the expected ro-independence of the mode when properly normalized (see 
([77])). In the following table we show the first eigenvalues corresponding to even boundary conditions at the 
tip. 







n= 1 


3.450 


n = 3 


22.113 


n — 5 


52.325 


n = 7 


94.558 


n = 9 


148.845 



In section [5] we prove the stability of the configuration by transforming the differential equation (|76| into a 
Schrodinger like one (see appendix [A]) . 



4.2 Non-Extremal D3-branes 

The in-plane fluctuation equation of motion lj22j) for the background (|43]) takes the form 
(P 4 -P 4 f* d 



Vp 4 -1 dp 



(/-l)i 



Sxi(p) = 0, 1< j oq<P<oo. 



(80) 



The background i!43tl asymptotes AdS and therefore the asymptotic (p — * oo) behavior of the solutions of 
([801 is given by ([79]) . As in the case of the last subsection, we expect to find a discrete set of eigenvalues 
leading to normalizable solutions. 

We plot in figure[2Q]the lowest eigenvalue we found when numerically solving |80|) looking for normalizabe 
solutions. A zero mode appears precisely at the critical point of the length function L(po), that is for 
po ~ 1.177 (see fig. [4]), in agreement with l[75j) . 

We conclude that the (left) branch in figure [4] having L'(p ) > is unstable under linear perturbations. 
Finally, note that the numerical analysis indicates that solutions on the (right) branch in figure [4] having 
L'(po) < are stable under linear perturbations. Nevertheless, as discussed at the end of section 1531 one 
expects the solutions with 1.177 < po < 1-524 to be metastable, decaying to a pair of free quarks. 
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Figure 20: Lowest uo 2 eigenvalue of ([80|l giving a normalizable solution as a function of po. A zero mode 
appears for po ~ 1.177. The classical solutions with po < 1.177 are unstable against linear perturbations. 



4.3 Maldacena-Nunez background 

The in-plane fluctuation equation of motion for the Maldacena-Nunez background lf47j) takes the following 
form (El 



d_ 

dr 



'( e 20M _ e 20(r o ))f d 



o20(r) 



dr 



5x\ (r) = , < ro < r < oo . 



(81) 



where uu 2 = uo 2 a' N. We now compute the asymptotic behavior of (f83~|) to see whether we should expect a 
quantized spectrum or not. In the r — > oo limit the equation of motion for fei(r) reads 



dr 



i d 



— e r 4 — \ + ui e r 4 



dr 



-.2 „r„ 



r > 1 



(82) 



where we used that e 2< ^ r ) — * e 2l Y 2 for r 3> 1. This last equation can be written as 

dr 2 



4r ) dr 



(83) 



Sxi(r) = 0, r > 1. 

The r _1 term can be omitted in the large r limit and the asymptotic solution to l[81j) is then 

fei(r) ~e-^(/3 e ra + /3 ie - TO ), r » 1 , 

where a = V1 ~ 4 "" . Fr om (|83|) it follows that any > lead to normalizable solutions, the spectrum 
of stable in-plane fluctuation is therefore continuum. In the Cu 2 < case (a > |) Po must be set to zero 
and we have the possibility of getting a discrete spectrum of negative eigenvalues. Our numerical analysis 
could not find any normalizable negative eigenmodes, suggesting the stability of the classical configuration 
in agreement with the concavity condition lfl7|) . 

In the section[5]we show that negative eigenvalues does not exist from the study of a Schrodinger equation 
analysis of (|8lj) (see appendix |A|) . 

4.4 Klebanov-Strassler background 

The equation of motion for the in-plane fluctuation in this case takes the form 



d ( K(r) 
dr \ h(r) 



Mr) 

h(ro) ) dr 



6K(r) 



h(r) 
h(r ) 



Sxi (r) = , < tq < r < 00 . 



(84) 
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with to 2 dimensionless and K(r) and h(r) given by ([561 and f59|) respectively. From the r — ► oo limit of 
if (r) and /i(r) one obtains 



which gives 



dr 2 



d /e r d 
dr \ r dr 



1 



2l 



dr 



Jxi (r) = 0, r > 1 



5xi(r) = 0, r > 1 . 



2^ 



(85) 



(86) 



In the large r limit the r 1 and the last term in ([86]) can be omitted and the in-plane fluctuation asymptotics 
turns to be 

6xi (r) ~ a + aie" r , r > 1 . (87) 

The integration constant ao must be set to zero to obtain normalizable solutions and we therefore expect 
to get a discrete eigenvalue spectrum. In figure lf2"T|) we plot the lowest eigenvalue of |84|) we have found 



Figure 21: Numerical solution for the lowest w 2 leading to a normalizable solution as a function of ro for the 
Klebanov-Strassler background. No negative eigenvalues where found. 

numerically leading to a normalizable solution as a function of ro. We have not found numerically any 
negative eigenvalues. In section [5] we will prove the stability of the classical solution transforming lj84|) to a 
Schrodinger like equation and showing that no negative modes can appear (see appendix [A} . 



4.5 Generalized Maldacena-Nunez 

The equation of motion for the in-plane fluctuation Sx\(p) in the backgrounds f63|) takes the form 



d_ /( e 8/(P)- e 8/(Po))| d 
dp [ 



D 8/(p) + fc(p) 



dp 



Sxi (p) — , < po < p < oo . 



In the p 



case (k(p) — const.) the equation f81j) for the Maldacena-Nunez background is recovered 



(from now on we consider fj, ^ — |). In the large p limit the gMN solutions asymptote the deformed conifold 
and the / function approaches a constant /oo, the asymptotic behavior is then given by 



-Hp). 



-Hp). 



Mo. 



dp \" dp J ' ~ e 8 ^~ — e s f( po 
Returning to the original r variable in |63|l (dr = e k ^dp) one obtains 

Td 2 



5x 1 (p) = 0, p>l. 



dr 2 



■Q 2 



5xi(r) = 0, r > 1 



(89) 



(90) 
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Figure 22: Lowest u 2 eigenvalue of ([88]) lead- 
ing to a normalizable solution as a function of 
Po- Negative (unstable) modes are found pre- 
cisely for the classical embeddings satisfying 
L'(po) > 0. Colors as in fig. M 



Figure 23: Zoom of figure [22] near the ori- 
gin. The lowest eigenvalues are positive for 
L'(p ) < solutions. A zero mode appears 
precisely for the critical values of the length 
function L(po) (see fig. [T4]) . 



whose solutions are plane waves e ±lu)r for cD 2 > and real exponentials e ±ujr for Co 2 < case. We conclude 
that no normalizable solutions exist for lj 2 > 0. A word of caution, as discussed in section [375l the gauge theory 
brane must be placed at a finite distance p^ , therefore for f88|) defined on po < p < p^ positive eigenvalues 
will exist. In the uj 2 < case the possibility for negative eigenmodes exists and in fact we find normalizable 
negative mode solutions precisely for the classical solutions region where the convexity condition (fl7| is not 
satisfied. In the figures l22l and l23l we plot the minimal eigenvalues leading to normalizable solutions we found 
numerically as a function of r . We found complete agreement with figure [HI no instabilities are found for 
classical solutions satisfying L'(ro) < 0, on the other hand, we find negative (unstable) modes for the right 
branch curve (Z/(ro) > solutions) in figure [14| These results are gratifying since unstable modes are found 
precisely for the classical embeddings which do not satisfy the conditions lfl7|). In the following section we 
review this results by transforming the equation into a Schrodinger like problem. 



5 Schrodinger Potentials Analysis 

In this section we analyze the fluctuation equation of motion (|22|) transforming it to a Schrodinger like 
equation (see appendix[Aj . From the form of the potential it is possible in some cases to show that no negative 
eigenvalues can appear and therefore to prove the stability of the corresponding classical embeddings. 



5.1 AdS 5 

The Schrodinger potential (|107|) for the equation (|77| takes form [15] 



V(p) 



P 4 -l 
P 2 



p e [i,oo) 



(91) 



here p should be understood as p = p(y). The change of variables (| 105^ leading to the Schrodinger equation 
(|106p can be analytically computed 



y(p) = Vo 



(92) 



with y = iW=. The half line p G [1, oo) of the original Sturm-Liouville problem, under the change of 
variables lf92j) . maps to the finite interval y € [0, yo], the potential lf9ljl diverging at yo- We have therefore 



obtained a Schrodinger problem defined on a finite interval with canonical boundary conditions (see (| 108[i - 
(|109p ) hence a discrete spectrum will result, moreover, since the potential (|9"Tj) is positive definite a standard 
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Figure 24: Schrodinger potential (|9Tj) for in- 
plane fluctuations in AdS the positive defi- 
nite property of it guarantees that no negative 
eigenmodes exist. 
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Figure 25: V(p, po) (eqn. l|93|0 for different po 
values. The blue line corresponds to po = 1.1, 
the green to the critical value po = 1.177 and 
the red, black lines to po = 1, 2. The region 
where the potential is negative diminishes as 
po increases. Negative eigenmodes cease to 
exist for po > 1.177 



QM argument tell us that no negative eigenvalue solutions exist. We conclude that the AdS embedding 
given by lj34j) is stable under linear perturbations. Figure [24] shows a plot of the potential l|9ip as a function 
of p, the true variable for the Schrodinger problem is y given by ([92)1 and it amounts to a rescaling of the 
horizontal axis in fig. [24] mapping p = oo to a finite distance. 

5.2 AdS'5-Schwarzschild 

The Schrodinger potential (|107|) for the in-plane fluctuation equation l|80p takes form |15j 

,P 8 (p 4 -Pe!)-^(4p 4 -l)-3p 4 



V(j>,Po) 



P 6 (P 4 1) 



1 < Po < P < °o 



(93) 



The behavior of this potential for different values of ro is shown in figure [25] Unlike the AdS$ x S 5 case, 
there exist regions where the potential becomes negative, this is in agreement with the results of section 
4.21 where negative eigenmodes where found. The potential starts from a negative value at po given by 
V(po,po) = —8/pq. As po increases the negative region gets dimmer and the negative modes cease to 
exist at the critical value, found numerically in section |4~2] p 0c — 1.177 which precisely coincides with the 
critical value of the length function L(p ). We conclude that the classical embeddings satisfying L'(p ) > 
are unstable under linear perturbations (see also [U] for a perturbative analysis of the eigenvalues). The 
classical unstable solutions L'(po) > have regularized energy E q q greater than zero (see fig. [6]), since the 
reference configuration satisfies the same boundary conditions, the natural candidate for the decay process 
is the reference (free quarks) state. 

For completeness we quote that since the asymptotics of this background coincides with the previous 
case, the Schrodinger equation for in-plane fluctuations results defined on a finite interval. The spectrum of 
stable fluctuations is therefore discrete. 



5.3 Maldacena-Nunez 

The Schrodinger potential for lj81j) takes the form 

-200) 



V(r, r ) = "—^ ([e 2 ^ r) - 3e 2< ^° VV) + 2(e 20W + e 2< ^° V(r)) , < r 



< r < oo . 



(94) 
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Figure 26: Schrodinger potential i|94|) . The 
blue, black, red and green lines corresponds 
to ro = 0.2, 0.7, 1, 2. The minimum of the 
potential becomes negative for ro > 1.1605. 



V(r,r„) 




Figure 27: Potential (|94|) as a function of ro. 
The potential asymptotically tends to Vx = j 
in agreement with l|83p . For ro < 1.1605 is 
positive negative and for r > 1.1605 contains 
negative regions. Near the tip the potential 
can be approximated asymptotically (for large 
ro) by V(r) ~-± + §(r-r ) 



As before r should be understood as r = r{y) and contrary to the last two cases the change of variables (| 105[) 
give a Schrodinger problem in y coordinates defined on the half line y € [0, oo). Figures l26l and l27l show 
the Schrodinger potential |94|) for different values of ro. We should confront these figures with the results 
in section l4~3l The outcome of that section, for all values of ro, was that a continuum spectrum results 
for u> 2 > and numerically no negative normalizable modes were found. We first address the continuum 
spectrum for lj 2 > 0. Figure [26] shows that the potential is positive definite for ro < 1.1605 and asymptotes 
the value Voo — j- One might therefore conclude no solutions for < Hi 2 < V m in, a discrete spectrum for 
V m in < u) 2 < \ (if possible) and a continuum, but not normalizable, spectrum for uj 2 > j, all in contradiction 
with the mentioned results. The agreement is achieved when taking into account the factor (PQ)~i that 
relates the solution of the Schrodinger equation with the fluctuation 5x\ (see appendix [A] eqn. I|105p ) 

5xi = ; ; r*~e _ 5*, r -> oo . (95) 

( e 20(r) _ e 2(Hr ))! V ' 

The e~5 factor makes all oj 2 > solutions of the Schrodinger problem satisfy the 5xi\ r =oo = whether or 
not they normalizable as ^(y) solutions (asymptotically one has y ~ r). However, for oj 2 < solutions the 
factor is not enough for making the (diverging) solutions satisfy the boundary condition. We conclude that 
for all ro a continuum spectrum results for Co 2 > 0. 

The remaining point to be addressed is the possibility of bound states for — \ < uj 2 < in the limit 
of large r . As seen from figure l27l asymptotically, the potential starts from V(r ) — —\ and the linear 
approximation one obtains is V(r) ~ — \ + |(r — ro). The relation between the r and y coordinates (|105p 
in the same limit is (r — ro) ~ y 2 /2. All these leads to a harmonic oscillator in y coordinates with bound 
state energy above zero. We therefore conclude that no bound states exist. We find a complete agreement 
between the Schrodinger analysis and the numerical results of section 14.31 

The results are appealing since if instabilities were found, no obvious candidate for the decay is available 
(cf. section l5~5j) . 
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Figure 28: Schrodinger potential ([96| . The blue, black, red and green lines corresponds to r = 0.2, 0.7, 1, 2. 
The potential is positive definite and therefore no negative (unstable) eigenvalue modes result. 
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Figure 29: Schrodinger potential (|97j) for dif- 
ferent po values and fi = — 1. The blue, 
black, red and green lines correspond to po = 
0.2, 0.7, 1, 2 respectively. 



Figure 30: Schrodinger potential for po = 0.2 
and different values of the parameter p. The 
blue, red, brown and green lines correspond to 
H = -1.8, -1.5, -1, -0.8. 



5.4 Klebanov-Strassler 

The potential for the Klebanov-Strassler in-plane fluctuation (|84| reads 

nr.ro) = - 8h3{r)h(ro) [4h(r)(h(r)+h(r ))ti(r)k'(r) 

-k(r)[3h{r) + 7h(r ))h' 2 (r) + 4h(r)(h(r) + h(r ))h" (r))] (96) 

The asymptotic behavior of the P, Q functions in this case lead to a Schrodinger problem defined on a finite 
y interval (see appendix [Aj . Figure |28|) shows the form of the potential for various ro values. The potential 
is positive definite and therefore no (unstable) uj 2 < solutions exist. The finite interval on which the 
Schrodinger problem is defined implies a discrete set of eigenvalues. We find complete agreement with the 
results of section 14.41 

5.5 Generalized Maldacena-Nunez 

The Schrodinger potential for in-plane fluctuations in the generalized MN solutions 1|88|1 can be written as 

V ^ f*) = e8 /(p) 2 +2fc (p) ( 2 ( e8/(P) " e8/(P0) )/' 2 W + (e 8/(p) + e 8/(po) )(/"(p) - k'(p)f(p)j) (97) 
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The asymptotic behavior of the P, Q functions (see appendix [A] and eqn. (J89J) ) leads to a Schrodinger 
problem formulated on the half line y G [0, oo). Figure (|29| shows the behavior of the potential for different 
values of po and a fixed value of p = —1. The potential becomes negative above some critical value p» 
and asymptotes V^o — in concordance with lj89j) and the existence of negative (unstable) modes found 
numerically in section l4~5l The (PQ)~* factor relating the Schrodinger wave function ^ to the fluctuation 
Sx i approaches a constant at infinity, therefore not changing the asymptotics of the 4" solutions (cf sect. 
15. 3j) , Figure [30] shows the Schrodinger potential, for a fixed value of po = 0.2, for different values of p. The 
minimum of the potential decreases as the p approaches — |. As already mentioned the original MN solution 
((47]) is not continuosly connected with the generalized class of solutions lj63|) . Agreement with the numerical 
results of section 15.51 is found but it is no clear to us which is the final state of the decay. 

6 't Hooft loop 

The electromagnetic dual to Wilson lines in Yang-Mills theories are the 't Hooft lines [35]. In four dimensions, 
the mechanism for confinement is supposed to be due to magnetic monopole condensation (dual Meissner 
effect), the analysis in [35j concluded that a screened monopole potential between a mm pair should be 
observed when confinement is due to a dual Meissner effect. A generalization of this idea is dyon confinement 
and goes under the name of oblique confinement. 

The string prescription for computing 't Hooft loops in the MN and KS solutions was proposed in the 
same papers [25|-[26] (see [36] and also [31] where a technical issue, correcting the proposed 2-cycle in [25] , 
was pointed out) and consists in wrapping a probe £>3-brane on the same 2-cycle on which the £>5-branes 
leading to the backreacted geometry were wrapped (see also |9j). The outcome of the construction is an 
effective Dl-brane (string) which is analyzed in complete analogy with the probe fundamental string we have 
been discussing in previous sections. The important difference with respect to the Wilson loop case is that 
the 't Hooft loop in generic non S-dual theories is sensible to the internal five dimensional manifold. 

6.1 Maldacena-Nunez 




Figure 31: The length function for the ef- Figure 32: The energy of the monopole- 

fective string as a function of r in the antimonopole pair as a function of the sep- 

Maldacena-Nunez t Hooft loop case. aration length. 

The embedding manifold for the D3 in the metric f47j) is |36p°l 

Ma = [t, x, r(x),9 = 0, ip = 2tt - <p, ip = tt] . (98) 

10 The remaining coordinates are set to constants. The value of the ip coordinate is fixed by demanding the S 2 to be of 
minimal volume. 



2fi 



The induced metric on VW4 results in 



'ind 




(99) 



and the expressions f49|) for a, h give 




(100) 




(101) 



Seff = 4nT D3 T(a'N) 2 / r tanhrVl + r 2 dx . 



(102) 



The important difference wrt the previous Wilson loop calculation in the Maldacena-Nunez background (see 
sect. I3.3P resides in the f(r) and g(r) functions in l|50p - lj5Tj) being multiplied by the 2-sphere volume (|100|) . 

In figures [31] we plot the behavior of the length function |(9|) as a function of ro . The length function 
is an increasing function of ro and from the previous discussions we therefore expect the embedding to be 
unstable. The instability of the embedding can be easily seen in the Maldacena-Nunez case since a fluctuation 
along the si-direction depending only on t, r is decoupled from the angular ones (consistent fluctuation). 
The Sxi equation of motion results in ((22]) with f(r) — g(r) — h(r) — r tanhr e^ r \ The asymptotics of 
the fluctuation is the same as in the Wilson loop case, nevertheless the behavior drastically changes near 
the origin since /(r) goes to zero. As seen in figure [33] negative eigenvalues exist for all ro values. For 
completeness we plot in figure l34l the Schrodinger potential associated with the in-plane fluctuation equation 
of motion. 

In figure [32] we plot the energy as a function of the endpoints separation length L. The energy of the 
configuration is positive for all L, this fact and the instability of the configuration suggests that the stable 
configuration for given boundary conditions is the one corresponding to two "straight lines". Contrary to 
the Wilson loop case the "straight lines" (used as reference state for regularizing the energy) can end at the 
origin since they correspond to wrapped D3 on the topological S 2 of {47]) which smoothly collapses at the 
origin. 

6.2 Klebanov-Strassler 

This case again involves wrapping a D3 over a topological S 2 inside f54|) (for its parametrization see the 
appendix A of [37J). There are important differences with respect to the MN case, in the present case the if 3 
supporting the KS geometry contributes to the string action (Q]) , and moreover it leads to the entanglement 
of the angular and the in-plane fluctuations. The different UV behavior wrt the MN is the reason for the 
length function L(r ) having a priori stable regimes (see figure 135]) . The behavior of V(L) shows the potential 
is screened for large L and this agrees with the linear confinement potential for the Wilson loop case (see [9] 
for a related example). The analysis of the present case is analogous to the one in sect. 13.21 summarizing 
when the configuration energy becomes positive, the two "straight lines" solution becomes favored. As for the 
MN case this becomes possible without a horizon in this case due to the D2> being wrapped over a smoothly 
collapsing S 2 at the origin. We did not attempt the analysis of the coupled fluctuations equations of motions 
to check for instabilities on the left branch of fig. l35l 

11 Placing the gauge fixed ansatz (198l> into the action (1 10 111 give the correct equation of motion for r(x) which coincides with 
0. 
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Figure 33: Lowest numerical eigenvalue J 1 for 
in-plane fluctuations depending only on t, r co- 
ordinates as a function of ro in the Maldacena- 
Nunez 't Hooft loop case. 




Figure 34: Schrodinger potential for in-plane 
fluctuations depending on t, r-coordinates as a 
function of ro in the Maldacena-Nunez t Hooft 
loop case. Blue, black and red lines corre- 
spond to ro = 0.2, 0.7, 1. 




Figure 35: Length function L(r ) for the ef- 
fective string as a function of ro in KS back- 
ground. 



Figure 36: The energy of the monopole- 
antimonopole pair as a function of the sep- 
aration length. 



7 Conclusions 

In this work we have analyzed the string proposal for computing rectangular Wilson loops via string embed- 
dings in gravity backgrounds and we have studied their stability under linear perturbations. 

The string prescription involves solving for a minimal open string worldsheet whose endpoints lie on the 
loop to be computed located at a fixed value of the holographic radial direction. When the endpoints are 
moved to infinity a divergent area results and a regularization is mandatory in order to get a meaningful 
answer. In section [2] we have reviewed this prescription and showed how a finite value is obtained. We have 
chosen to regularize the action by the standard procedure originally proposed in [4]. This is interpreted as 
saying that the Nambu-Goto worldsheet area computation includes the interaction energy plus the self energy 
(mass) of the external quarks. Within this interpretation we reproduced the well known results for AdS and 
thermal AdS. The regularization was in fact responsable for turning the original positive area into a negative 
attractive potential energy. When turning to smooth backgrounds (AdS in global coordinates, MN and KS) a 
puzzle arises since the straight strings running along the radial direction used in the substraction prescription 
must end somewhere in the bulk. We concluded that the correct interpretation for the substraction procedure 
is that we are comparing the string 'Wilson loop' worldsheet with respect to a reference state consisting in a 
straight string worldsheet whose endpoints lie at antipodes of a compact direction (pictorically represented 
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in fig. [TO]) . It then follows that the reference state in general satisfies different boundary conditions than the 
worldsheet used in computing the expectation value for the rectangular Wilson loop. This last observation is 
welcomed for the MN and KS cases where the linear confining relation occurs for worlsheets having positive 
regularized energies (see figs. [9] and [13]): if the reference state satisfied the same boundary conditions as 
the 'Wilson loop' worldsheet, the observed linear behavior should not be considered since the reference state 
(Egg = 0) would have been the lowest energy one (cf. last paragraph of sect. 13 . 2[> . but from the previous 
analysis we see that this is not the case. We would like to recall an observation in [9J stating that the relation 
between Wilson loops and strings in gravity duals (at the semiclassical level) 

(W) ^ e~ A (103) 

is schematic since the addition of boundary terms to the Nambu-Goto action does not change the minimal 
area character of the solutions but turns the value of the classical action into something different than the 
area. In [9] this arbitrariness was used to make a Legendre transform of the Nambu-Goto action showing 
that the resulting quantity, for the case of loops in AdS, is free from the linear divergences arising from the 
behavior of the worldsheet near the boundary of AdS. 

We also discussed the concavity conditions ifTTj) that must be satisfied by any potential pretending to 
describe the interaction between physical quarks. Generic gravity duals have positive and increasing f(r) 
functions, so the concavity conditions are not satisfied when the length function is an increasing function of 
the minimal radial position reached by the string r . In section [3] we performed the analysis of the length 
and potential functions L(tq) and V^ tr i ng (i) for different backgrounds and showed that some of them lead to 
embedding solutions where the concavity condition fails. 

Based on previous work [14]-[T5] we studied linear fluctuations around the embedding to test the stability 
of the classical embedding. We concluded that whenever the solution leads to an unphysical potential, not 
satisfying the conditions ifTTI there exist unstable modes under linear fluctuations. In the course of the 
analysis we discussed the different gauge fixings that can be imposed and its relation with the diffeomorphism 
of the Nambu-Goto action. Three natural gauge fixings where discussed and we chose to work in the r-gauge 
since it lead to simpler closed expressions for the fluctuations equations of motion (see eqn. l(2lj) ). The 
r-gauge leads to singular behavior in the fluctuation at the tip of the embedding, but reviewing [15] we 
showed that nevertheless they are physical once an appropriate gauge transformation is performed. 

In section |4] we perform the stability analysis for the solution reviewed in section [3l We showed by a 
numerical analysis that the AdS§ x S 5 , Maldacena-Nunez and Klebanov-Strassler are stable. On other hand 
for thermal AdS and the generalized Maldacena-Nunez backgrounds of sect. I3.5l we found unstable modes in 
agreement with the behavior of the L(ro) relation. This last case is rather pathological since the loop cannot 
be placed at infinity and moreover we found that a minimum separation exists beyond which no smooth 
solution connecting the string endpoints exists. In section [5] transforming the Sturm-Liouville fluctuation 
equations of motion into a Schrodinger like equation we reanalyzed the problem, finding complete agreement 
with the results obtained in section |4j We conclude that the regions where we find unstable modes coincide 
with the regions where the concavity condition fails. 

In the last section we performed the previous analysis for the case of monopole-antimonopole interaction 
in the non-conformal gravity duals of Maldacena-Nunez and Klebanov-Strassler. We discussed the 't Hooft 
loop string prescription given by wrapping a D3 on the topological S 2 present in the geometries. The MN 
case was shown to be unstable for all ro values. A fluctuation analysis was feasible since a decoupled equation 
for the in-plane fluctuation could be found were an unstable mode was shown to exist by a numerical analysis. 
The KS presented a behavior similar to thermal AdS with presumably stable and unstable regions, but a 
the fluctuation analysis lead to coupled fluctuations equations of motion which we did not analyze. 

We conclude that the analysis of Wilson/ 't Hooft loops in given gravity background by looking at the 
value of the / 2 at the origin should be supplemented with an analysis of the L(ro) relation. 
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A Sturm-Liouville to Schrodinger 

Equations l|22p - l(2"3|) are of the Sturm-Liouville type 



^- ( p ( r , r o)-^j- ) +U(r,r ) $(r) = w 2 Q(r. / (1 K>( /■ ). ,•„ ; / x 



(104) 



the functions P(r, ro) and Q(r, ro) can be read off from ([22"l) - l|2"3l) . U(r, ro) = in both cases. The change of 
variables 



y 



transforms <|104[) to a Schrodinger like equation 

d 2 



dy 2 



V 



(105) 



(106) 



Here yo — dry j; which may be finite or infinite depending on the nature of Q, P and one can check that 
(|105p is integrable at the lower limit giving y ~ \Jr — Fq. The potential V is given by 



V = — 



{PQY 



Pi d I IP d 
Qidr \VQd^ 



{PQY 



(107) 



The points r = ro and r = oo map to y — and y = yo respectively. The boundary conditions to be imposed 
on the solutions of (|106[) are |15| : 

• Infinity: string endpoints fixecf^l 

Sx\ r=oc = => y\ y=yo = (108) 

• Tip r = r : for both in-plane 8x\ and transverse fluctuations Sx m one obtains from (|24| . (|3"Tj) . i[32|) 

= 



Even solutions : 4^ 

dy 



Odd solutions : \&\ y=Q = . 



(109) 



B Exact spectrum for transverse fluctuations in AdSs x S 5 

We review here the solution of [13] for the exact spectrum of the longitudinal fluctuations in the AdS§ x S 5 
background and compare it with our numerical results using the shooting technique described at the end of 
section 12.21 

The AdS metric is written in Poincare coordinates 



K 

ds 2 = -^-{-dt 2 + dxidxi + dz 2 ) + R 2 dnj 



A x-gauge fixed ansatz t = r, x = x c \ , z — z c \ (x) leads to 



dz d 
dx 



fe,) 4 



(110) 



(in) 



12 See however a loophole in the Maldacena-Nufiez context (sect. 15. 3|l when imposing l|108 
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The solution to with the string endpoints separated by a distance L is ([34|) - (|35"1) . 



Xc\{z) = ±z 



2r[i] 2 4 V^o'4'2 



(112) 



where zq = z c i(0) = (r[^] 2 /(27r) 2 )L is the maximal radial distance reached by the string (tip of the string). 
Fluctuations around the solution (|1 12[) in the transverse x m (m = 2, 3) coordinates decouple, writing 
= (t 7 x c \(<j),Sx m {t,a),z c \(cr)) the equations to linear order are [13j.[T7] 



x— gauge : 



r— gauge 



q2 z d ( x ) q2 



5x m {t,x) = 



d*-{l--)dl + -d z 



Sx m (t, z) = m = 2, 3 . 



(113) 
(114) 



As mentioned in section l2Tl note that the x-gauge equation of motion (|113[) depends explicitly on the classical 
solution z c \(x). The equations are related by the change of variables given in (j 1 1 2j) . Writing 5x m = e~ %wt f{z) 
in l|114p and calling z = z/z one obtains [17] 



(i-m 2 -§ +e 



f(z) = 0, 0<z<l, 



(115) 



where £ = zqlo. The change of variables |38J 

f(i) = yf+WFM 



„l) = ±2 /'- 



(l + (^)2)^T3t4 

transforms equation (j 1 1 5 j) into a simple harmonic oscillator 

^f + I^ 4 - 1 ) F = > ?e[-9..9.1 



(116) 



(117) 



where g* = g(0). The boundary conditions at infinity <5:r m (i,0) = have been mapped to -F(g*) = 0, and 
quantize the frequencies in l|117p leading to 



U! n Z \ LU n Z Q 



Zt\ 1 



t 2 dt nw 
= — . n = l,2,. 

(i + w*4)VT=F 2 



(118) 



The following table shows the comparison between the exact eigenvalues (|118p and our numerical calculation 
of eigenvalues of (| 1 15|) with zq = 1. 





Exact 


Numeric 




2.203 


2.226 




3.467 


3.492 


W 3 


4.697 


4.735 




5.914 


5.959 


W 5 


7.125 


7.181 


W 6 


8.332 


8.396 




9.537 


9.612 


u>$ 


10.741 


10.823 



The odd (even) eigenvalues where obtained solving (|115|1 with the even (odd) boundary conditions discussed 
after (El. 
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